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ABSTRACT 

Using COBE-DMR 4-year data, we find a general expression yielding the likelihood dis- 
tribution in the 3-dimensional parameter space spanned by the spectral index n, the spectral 
amplitude a 10 and the false-vacuum density parameter Q^. Using such simple expression, the 
range of possible normalizations, within a given likelihood interval from top-likelihood normal- 
ization, is readily found, with fair approximation, for any model with total density parameter 
Q = 1 and assigned n and Q\. 



Subject headings: cosmic microwave background 
ture of the Universe - methods: numerical 

1. Introduction 

The COBE-DMR data, so far, provide the only 
all-sky map for the cosmic microwave background 
radiation (CMBR) (Smoot et al. 1992). Using 
such data, the angular temperature fluctuation 
spectrum Ci can be inspected for I values from 2 
to 30. In principle, cosmological model parameters 
should be adjusted so to provide the best possible 
fit to all Ci. Unfortunately, however, the situation 
is more complicated than so. As a matter of fact, 
model independent Ci values and error bars were 
provided by Hinshaw ct al. (1996), Bunn & White 
(1997), Tegmark & Hamilton (1997), assuming 
that 1(1 + l)Ci takes single values, in suitable I 
intervals. Tegmark (1997) and Gorski (1997) also 
provided angular spectrum data and error bars, 
separately for all I values. Attempting to use such 
Ci to fix model parameters, as outlined also by 
Gorski (1997), is far from trivial. In fact, two main 
difficulties arise. First of all, each C; distribution, 
around its maximum likelihood value, is not Gaus- 
sian, as is witnessed by its non-symmetric 1 a er- 
ror bars; deviations from a Gaussian behaviour are 
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even stronger beyond 1 a. Furthermore, the C\ es- 
timators are strongly correlated with each other. 
Therefore, even obtaining the likelihood distribu- 
tion for each C;, up to a few a' s around their max- 
imum likelihood values, a simultanous direct fit to 
different C\ 's bears no clear statistical significance. 
Then, if precise likelihood comparisons within a 
given class of models are to be made, one has lit- 
tle alternative to fitting directly model parameters 
to COBE-DMR pixel data. Such data are how- 
ever public and can be found on the NASA site 
(http://space.gsfc.nasa.gov/astro/cobe/). 

In this work we focused on models with total 
density parameter Q a = 1, due to contributions 
of baryons (density parameter Q&), cold-dark- 
matter (CDM; density parameter f2 c ), and false 
vacuum (density parameter Q\ — A/8irGp cr ). Re- 
cent outputs of BOOMERanG and Maxima I ex- 
periment (de Bernardis et al. 2000, Hanany et al. 
2000) confirm that f2 Q values substantially differ- 
ent from unity are disfavoured. fi c and to- 
gether with the primeval spectrum index n and 
the primeval spectrum normalization, are the pa- 
rameters which mostly shape the C; spectrum. On 
the contrary, up to I ~ 30, the dark-matter (DM) 
composition and the Hubble constant H bear a 
limited impact only. Attention was recently con- 
centrated on models with A / 0, mostly because 
data on SN la (Riess et al. 1998, Pcrlmutter et 



1 



al. 1999) seem to favour a negative deceleration 
parameter. 

Previous analyses of the likelihood distribution 
in the ai$-n plane were performed by various au- 
thors. The parameter aio = VCio Tq accounts for 
primeval spectrum normalization in a simple way, 
as its best-fitting value is nearly independent of n 
in the 4-year COBE data. 

Assuming a pure Sachs & Wolfe spectrum 
(Sachs & Wolfe 1967, Bond & Efstathiou 1987), 
Gorski et al. (1994), Bcnnct ct al. (1994), Wright 
ct al. (1994), Bond (1995) and Tcgmark & Bunn 

(1995) analysed the 2-year COBE data. Under 
the same assumption, Gorski et al. (1996, G96 
hereafter), Bennet et al. (1996), Wright et al. 

(1996) and Hinshaw et al. (1996), analysed the 
4-year COBE data. Such spectrum can be con- 
sidered a reasonable approximation for pure CDM 
or mixed models with fl Q = 1, for low Vs. On the 
contrary, open or A models are expected to behave 
differently (see, e.g., Stompor & Gorski 1994). A 
first attempt to determine the normalisation of 
A models and an upper limit on the density pa- 
rameter 17a using the 2-year COBE data was 
performed by Stompor, Gorski & Banday (1995a, 
b) and Bunn & Sugiyama (1995). But the most 
extended results on the likelihood of A-modcls, 
based on 4-year COBE data, were obtained by 
Bunn & White (1997, BW hereafter). 

Let us briefly summarize BW results, which 
were also partially included in the popular CMB- 
FAST code (Seljak & Zaldarriaga 1996), to show 
why a further effort is needed to complete them. 
BW define a function D(x) = 1(1 + l)C u with 
x = log 10 /, so that D(l) = 110 Cio, and consider 
the parameters D 1 and D" in the expansion 

D(x)^D 1 [l + D'(x-l) + ^-(x-lf]. (1) 

Accordingly, D' and D" essentially account for 
the first and second derivative of D(x) at x = 1 
(I = 10). BW provide, also through an analytic 
fitting expression, the maximum-likelihood nor- 
malization D\ for models defined by the values of 
D' , D" parameters and the likelihood distribution 
in D'-D" plane, for the top likelihood value of D\. 
Among 2-paramctcr fittings of COBE-DMR data, 
the BW approach is quite effective in treating A- 
models, as fixing D' and D" corresponds to giving 
n and A and viceversa. 



However, BW provide only a partial infoma- 
tion on the likelihood of models whose amplitude 
is close to (but not coincident with) the best-fit 
amplitude for given n. In fact, let n-aio yield the 
top-likelihood model and n± — n±An limit the 
1 a error bar along the n axis, as given by BW. 
They complete the information on likelihood dis- 
tribution, by stating that la error bars (on am- 
plitude) roughly correspond to 7%. 

In order to determine the likelihood of a model 
defined by n-a w values close to n-a w , one has 
to assume that, al least within 1 a, the likelihood 
distribution is Gaussian. If so, one may determine 
first the likelihood of the intermedate model given 
by n-aio. Then, still assuming a Gaussian dis- 
tribution along the aio axis, one may determine 
the further likelihood decrease when passing from 
n a 10 to n--a w . As a matter of fact, both distri- 
butions, along the n and the a w axes, are signif- 
icantly non-Gaussian, as is shown by the shape 
itself of the top-likelihood curve. This already af- 
fects the above procedure within 1 er, but makes it 
seriously unadequate above 1 a. 

It may be important, instead, to know how the 
likelihood is distributed among models with dif- 
ferent normalizations. For instance, varying the 
normalization of a model with given n and A, we 
may predict different galaxy cluster number den- 
sities. A model under-(over-)producing clusters, 
for maximum-likelihood normalization, might be 
in agreement with observations if a different nor- 
malization is taken. If this occurs for spectral in- 
dex and normalization within 1 a from the top- 
likelihood values n-aio, rejecting such model be- 
cause of cluster predictions is illegitimate. But 
also knowing the range of normalizations allowed, 
within 1-2-3 it's, for a model with given A and 
n may be relevant, as well as determining, for a 
given n and normalization, which is the range of 
A values allowed at the 1-2-3 a level. Similar ar- 
guments can be made for the capacity of a model 
to predict a fair amount of high-z objects or to 
fit the spectral slope parameter T, yielding the ra- 
tio between mass variances at 8 and 25 /i _1 Mpc 
(see, e.g., Bonometto & Pierpaoli 1998, for more 
details), while other parameter combinations may 
also have to be explored. 

Quite in general, to extend BW results in this 
way, we must perform a 3-parameter analysis of 
COBE- DMR data. Our parameters are the spec- 
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tral index n, the normalization aio, and the vac- 
uum density parameter Oa- The result of our anal- 
ysis will be a handable analytic expressions of the 
likelihood distribution in such 3-parameter space. 
We shall also provide a simple numerical routine, 
which can be appended to CMBFAST and gives, 
besides of the best-fit normalization for fixed n 
and A, also the limits of the fluctuation amplitude 
intervals where likelihood decrements correspond- 
ing to n ct's (n = 1,2,3) occur. Other details 
provided by such routine will be discussed below. 

2. From COBE data to the likelihood dis- 
tribution 

In this section we shall report how the likeli- 
hood distribution is obtained from COBE-DMR 
data. A large part of this section is based on 
Gorski (1994, G94 hereafter). Its details are 
unessential to read the following sections and to 
use our relations yielding the likelihood distribu- 
tion. Let us rather draw the reader's attention 
on the Appendix B, where we show in detail how 
the likelihood dependence on monopole and dipole 
terms are integrated out. 

COBE-DMR temperature data were estimated 
at the frequencies of 31.5, 53, and 90 GHz. Us- 
ing such estimates, the vectors d, with N p (pixel 
number) components di, can be built. The sig- 
nal variance, due to uncertainties of instrumental 
origin at each pixel, Ci, can also be collected in 
Afp-dimcnsional noise vectors (<x). Each pixel is 
centered on a point of the celestial sphere of co- 
ordinates hi = (6i,<j>i). The likelihood of a given 
model M. is obtained by comparing d with a fic- 
titious data vector S (components Si), built from 
the model, but taking into account how real data 
are obtained and, therefore, also the noise vector 
er. Noise correlation among pixels, expected to be 
small, will be neglected in our theoretical develop- 
ments. 

In principle this can be done by evaluating the 
matrix M, whose components 

Mij = (SiSj) (2) 

are obtained by averaging on the ensemble of M. 
realizations. Assuming a Gaussian statistics, the 
likelihood of M reads then: 

C = [(2^) Ar "detM]- 1 / 2 exp[-(l/2)d T M- 1 d] . 

(3) 



A slight non-Gaussian behaviour in data (Fer- 
reira, Maguejo & Gorski 1998, Novikov, Feldman 
& Shandarin 1998, Pando, Valls-Gabaud & Fang 
1998) possibly originated by post-recombination 
processes, should not affect our analysis, as shown 
in detail by Contaldi et al. (1999), for a standard 
CDM model. 

However, eq. (3) can hardly be used, because 
of the high dimensionality of the matrix M, which 
should be inverted. Replacing the pixel basis (a 
discretized coordinate representation) by an an- 
gular harmonic basis (essentially its Fourier trans- 
form) allows to reduce the dimensionality of matri- 
ces, without wasting physical information. When 
doing so, we must also take into account that tem- 
perature fluctuation data about the galactic equa- 
torial plane do not give cosmological information. 

In order to evaluate the likelihood of a model 
M, therefore, the first step amounts to defining 
the set of real spherical harmonics: 

Y lm (f) = VTTTJl^ $^j^ |m| (cos 9)f m (<f>) 

(4) 

where 

f m {4>) = (27T)- 1 / 2 , k- 1 ' 2 cos(m^) , 7T" 1 / 2 Bin (M0) 

(5) 

for m = 0, > 0, < 0, respectively, while the 

PT(x) = (-l) m (l - x 2 ) m ' 2 {d/dx)\ m \Pi{x) (6) 

are obtained by differentiating the ordinary Legen- 
dre polynomials Pi. Such spherical harmonics are 
built so to fulfill the orthonormality conditions 

/ d 2 rY lrn (r)Y Vm ,(f) = 5u,8 mm , (7) 

J Alt 

when the integration is extended on the whole sky. 
If this integration is replaced by a sum on the N p 
pixel centers, orthonormality is recovered by re- 
placing Yi m by 

Y lm (f) - wf x Y lm (f) (8) 

(a table of wf lx and directions on their computa- 
tion are given in Appendix A). 

Using such spherical harmonics, then, we can 
build the components of a fictitious signal vector, 
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for the point f i; yielding the i-th pixel, as follows: 



CMB OMR , noisc\V> \ 



1=0 m=- 



(9) 

Here, besides of afJ^ B , that are to be obtained for 
a realization of .A/f , we need the a"^ ISC originated 
by the instrumental noise and the window func- 
tion Wi >MIi for the observational apparatus (given, 
e.g., by Wright et al. 1994). Accordingly, as above 
outlined, the fictitious signal vector that we shall 
build will convolve the features of a given cosmol- 
ogy with the characteristics of the COBE-DMR 
apparatus. 

It is also convenient to put together the indices 
I and to, by defining 



(j, = l(l + l) + (m+l) . 



(10) 



If I values up to l max — 30 are considered, [i at- 
tains the maximum value \i m ax = {lmax + l) 2 = 
961. Accordingly, the functions Y^ define the pas- 
sage from the TVp-dimcnsional pixel basis to a 961- 
dimcnsional basis. In the sequel, latin indices shall 
be used to indicate vector components on the pixel 
basis only; accordingly, eq. (9) can be unambigu- 
ously rewritten as 



E( 



(11) 



(~ is used to indicate that window function effects 
are already taken into account). 

Let us then remind that the a™ B = a£ MB , 
that we shall build, depend on the model M and 
on its realization; ensemble averaging them, one 
would obtain the angular spectrum components 

n — l\r,CMB\2\ 
C l = \\ a lm I )• 

The use of spherical harmonics, however, gives 
place to a problem. When using real data, pixels 
within 20° from the galactic plane are to be ex- 
cluded (in the sequel, residual pixels will be said to 
belong to the cut-sky). First of all, this reduces N p 
from 6144 to 4016. Furthermore, in the cut sky, 
the set of polynomials is no longer orthonor- 
mal. Orthonormality, however, can be recovered 
by replicing the basis Y= {Y^} by the basis 

* = r y 

where T is a 961x961 matrix that we shall now 
define. In order to have (\E'\E' T ) C ., ! . = I ((....) c . s . is 



a product obtained summing on the set of pixels 
belonging to the cut-sky; I is the unit matrix), it 
must be 

r • (yy t ) c . s . ■ r T = i (12) 

and this is obtained if the matrix (YY T ) C S . is ex- 
pressed as a product L • L T and 



r = L 1 . 



(13) 



Here we shall use the upper-triangular matrix T, 
obtained by performing the so-called Kolewski de- 
composition, as done by G94. This choice is not 
unique; alternative possibilities are discussed by 
Tcgmark (1997). However, as we shall see, this 
choice allows a simple integration of the likelihood 
on monopole and dipolc data. 

The fictitious data vector 8, can be expanded 
on both Y and \& bases; let a M and c p be its com- 
ponents, respectively. Clearly: 



^ ] r * nvC v and — S ' L^ v c 



(14) 



Therefore, once Kolewski decomposition provides 
the L matrix, we readily obtain c M = (d\P M ) c . s .; in 
particular, besides of the CMB component c^ MB , 
the noise component c^ OIse is soon obtainable. Let 
us also outline that, thanks to the triangular form 
of the L matrix, the c M coefficients are linear com- 
binations of a v coefficients with v > fi. 
Using the basis ^f^, eq. (11) reads: 



^ = E c a<iv = E( c m 



CMB 



cf se )^ (15) 



and 



CMB 



Tij-noisc 



(16) 

is the new correlation matrix (assuming that signal 
and noise are both distributed in a Gaussian way 
and uncorrelated). Using M™ B and M™ iso , the 
expression for the likelihood of the model reads: 



C = [(2Tr) N det(M CMB + M. nmse )]- 1 ' 2 x 
x cxp[-(l/2)c T (M c ' MB + M noise ) _1 c] . (17) 

Here c is a vector with 961 components c M . 
Thanks to the triangular form of L, we can easily 
integrate the likelihood results on for fi = 1, .., 4 
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(such integration can be made in an exact way; 
details are given in Appendix B). Hence, in the 
cut sky, the effective basis dimensionality will be 
957 (for each signal frequency to be considered). 

To build the likelihood eq. (17), Ci are obtained 
using the CMBFAST program, convolved with the 
COBE-DMR window function. Then, according 
to eq. (14), 

M CMB = L T, {k CMB k CMB Ty L = jT^CMB^. 

(18) 

here C CMB is a diagonal matrix, whose diagonal 
terms are the components of the angular spectrum 
Ci, however is repeated 21+1 times (for all /x values 
which correspond to a given I). M. CMB is then to 
be set together with the noise matrix 

M noisc = ( c „oi S c c noisc T) = ft^^.^T ^ 

(the index , runs on the 4016 pixels of the cut- 
sky). 

In this way, the likelihood of a given model M. 
is finally evaluated and different models may be 
tested against COBE-DMR data. 

3. Results 

Using the approach described in the previous 
section, substantially coincident with G94 treat- 
ment, we considered a lattice of models for h vary- 
ing from 0.5 to 0.7 and ft a from to 0.8. For each 
model, the likelihood was evaluated for different 
values of n and aio- We used 4-year COBE-DMR 
data, performing a weighted average of 53 and 90 
GHz outputs; we verified that using the noisy data 
at 31.5 GHz does not add any substantial improve- 
ment. More precisely, in accordance with Tegmark 
& Bunn (1995) and Bunn & White (1997), we may 
set 

= ^(53)^" (53) + ^(90)^(90) 
ar 2 (53) + ^(90) 

and 

ar^ar\53) + ar 2 (90) (21) 

as data and noise vector components, respec- 
tively. Using such inputs, we worked out the 
3-dimcnsional curves C(n,a\o) for each model. 
As an example, in Fig. 1 we show it for a model 
with h = 0.5 and A = 0. From such curve, the 
isoprobability contours can be readily obtained. 



They encompass volumes corresponding to 68.3 %, 
95.4 %, 99.7% of the total volume below the 2- 
dimcnsional curve, respectively, and will be called 

1- 2-3 a contours, in the sequel. 

We found a fairly regular behaviour of ln(£) 
which can be expressed through a third degree 
polynomial in x= n — n° and y = a 10 — a° ; here 
n° and a° yield the peak position. More in detail, 
our fitting formula reads: 

-2 ln£' = Fix 2 + P 2 xy + 5.9j/ 2 + P 3 x 2 y + 

+P 4 ?y 3 + P 5 xy 2 - 0.7Sy 3 . (22) 

In Table 1 we give the values of the best-fit val- 
ues of the parameters (i — 1, ..,5), n° and aj 
for the models of the lattice. The expression (22), 
together with the values in Table 1, outlines that 
substantial deviation from a Gaussian behaviour 
in the error distribution are present; however, to 
put them under control, up to the 3 a level, it 
is sufficient to use third degree polinomials. The 
quality of the fits obtained in this way can be ap- 
preciated from Fig. 2, where we show how the 1- 

2- 3 a contours obtained from model analysis and 
fitting formula agree, for a set of typical cases. 
Residual discrepancies, visible for the 3 a contours, 
correspond to overall likelihood shifts ~ 10~ 6 . 

The fact that the expected non-Gaussian be- 
haviour is so simply fitted, is not the only finding 
of this work. In fact, we also find that all the coef- 
ficients Pi, n° and a° in the expression (22), can 
be fairly approximated using a single interpolating 
expression, as simple as 

a + bSlf + cflf . (23) 

In Table 3 we give the values of the interpolating 
coefficients a, 6, c and S. In Fig. B, we show how 
our fitting formula meets the Pi values (in two 
typical cases) , n° and a° obtained for the various 
models. 

In order to validate our algorithm, we pro- 
duced 120 CMBR sky realizations for cosmolog- 
ical models with assigned aio = 6.93 /j,K, Q\ = 
and n = 1.00, simulated their observation with 
COBE-DMR, and applied our algorithm to deter- 
mine the model likelihood in the space spanned by 
aio, f^A and n. The simulated CMBR sky was pro- 
duced using an algorithm based on the technique 
suggested by Muciaccia, Natali & Vittorio (1997). 
The Ci coefficients required were generated us- 
ing CMBFAST. The spectrum was then multiplied 
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h 


n A 


Pi 


P2 


Pz 


Pi 


P 5 


n° 




0.7 


0.8 


9.93 


1.14 


2.24 


-1.53 


1.72 


1.1711 


6.1579 




0.7 


12.31 


0.50 


2.36 


-0.71 


1.92 


1.1149 


6.2173 




0.6 


14.68 


-0.123 


2.17 


0.66 


2.07 


1.0831 


6.2547 




0.4 


17.46 


-1.014 


1.38 


3.28 


2.28 


1.0550 


6.2927 




0.2 


18.62 


-1.44 


0.89 


4.59 


2.37 


1.0489 


6.3073 


0.6 


0.7 


12.43 


0.42 


2.33 


-0.61 


1.92 


1.1138 


6.2176 




0.6 


14.69 


-0.168 


2.14 


0.73 


2.07 


1.0852 


6.2552 




0.4 


14.69 


-1.04 


1.36 


3.22 


2.28 


1.0597 


6.2944 




0.2 


14.69 


-1.46 


0.88 


4.49 


2.36 


1.0550 


6.3099 


0.5 


0.4 


17.41 


-1.052 


1.36 


3.27 


2.28 


1.0605 


6.2941 




0.3 


18.08 


-1.300 


1.06 


4.02 


2.33 


1.0574 


6.3039 




0.2 


18.45 


-1.458 


0.89 


4.47 


2.36 


1.0574 


6.3102 




0.1 


18.61 


-1.540 


0.77 


4.69 


2.38 


1.0598 


6.3136 




0.0 


18.33 


-1.58 


0.7 


4.5 


2.38 


1.0647 


6.3138 



Table 1: Values of parameter to fit the likelihood function, sec cq. (22). Such values are reported as an 
intermediate step of our work and mostly to show the residual small dependence on h of the fitting costants. 
Such dependence will be neglected in the final fitting expression. 



by the COBE-DMR experiment window function. 
COBE DMR features were used also to define pix- 
els and pixel noise variance. Using such fictitious 
data, we searched the top-likelihood point for all 
models in the n-a\o plane, using the same algo- 
rithm applied to COBE-DMR data. The effect of 
the 20° subtraction along the galactic plane was 
also tested, by searching the top-likelihood point 
for all models, also without galactic plane sub- 
traction. Our algorithm recovers n = 1.01 ± 0.13 
and aio = 6.92 [iK ±0.30 [iK (the latter standard 
deviation corresponds to ~ 4%). When the full 
sky is used, standard deviations are reduced by 
^20% and <~ 30%, respectively. For the sake 
of comparison, let us report that a more limited 
test, reported by BW, who made simulated maps 
with purely Sachs & Wolfe input spectra, gave 
An ~ 0.26 and a normalization discrepancy up 
to - 7 %. 

We also compared our results with those cases 
already treated in the literature. In Table 2 we 
present the results of previous works on n and Q, 
using 4-year COBE-DMR data. The last line of 
the table provides our estimate, assuming a pure 
Sachs & Wolfe spectrum. The dispersion in the 
values of the table is to be attributed to differ- 
ent choices of the "Galactic-plane cut" and differ- 
ent combinations of the 3 frequencies, while results 
should not be affected by different data compres- 



sion techniques. For each result, however, we in- 
dicate the method used by the authors. 

Here we shall report some details on the com- 
parison with the analysis of G96, on 4-year 
COBE-DMR data, assuming a pure Sachs & 
Wolfe angular spectrum, and with the outputs 
of BW. 

In G96, a (3*961) component signal vector was 
built setting together the data of the three maps at 
31.5, 53 and 90 GHz. An analogous (3*961) com- 
ponent noise vector is also built. Pixels only at 
galactic longitude above 20° are taken. Monopolc 
and dipole components of each map, which are not 
physically relevant to the power spectrum estima- 
tion, could be exactly removed by integrating over 
the first four components. We performed a similar 
analysis using (2*961) component vectors, based 
on maps at 53 and 90 GHz (the 31.5 GHz map, 
characterized by high noise level, was not taken 
into account). Results are shown in Fig. 4, where 
we report 1-2-3 a curves; the top-likelihood point 
we find is indicated by an empty triangle, while 
the G20+* top- likelihood point (see table 2) is 
indicated by a filled box. 

In BW, a weighted average of the signals at 
53 and 90 GHz is used, as we did in this work. 
The signal compression technique is however quite 
different from the one used here. The total number 
of pixels they use is 3890, after all pixels below the 
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n 




Q(j*K) 


Ql 




G20* 






1 


99 +0.24 

■ zz — 0.28 


16.71 + 3i2 


19. 


v — 1.25 








x 


9O+0.23 
ZO -0.29 


15 26 +3:93 


1 8 


04+1.25 

°^-1.20 


G30 + 


* 




1 


90+O.29 
— U.o4 


15.18+^6 


17. 


09+1.44 
— 1.34 


G20 + 


Z 


-g 


1 


11+0.38 
■ lx -0.42 


Ifi 3S+ 5:is 

lu,JJ -3 .69 


17. 


00 + 1.77 
•°°-1.68 


G30 + 


Z 


-g 





7Q +0.48 




16. 


c 7+1.92 
•°'-1.82 


G20 + 


* 


- 9 


1 


91 +0.24 
zl -0.28 


1 C 90 + 3.69 

10.ZO_ 2 64 


17. 


R7+1.25 
•°'-1.15 


G30 + 


* 


- g 


1 


9 4 +0.27 




17. 


J *-1.34 


H w 






1 


O r+0.26 

■ /o -0.29 


15.4^;g 


18. 


4 +1.4 
4 -1.3 


H io - 
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1 


90+O.26 
■ ZO -0.27 




17. 


o + 1.3 
•°-1.3 


H J 






1 


nf) +0.40 


17.2tS:o 


17. 


9+1.9 
• z -1.7 


B w2 






1 


.18 ±0.28 


16.2 


18. 


.7 ± 1.26 


T20 w2 




1 


91 +0.35 
• zl -0.42 


1 4 y+5.56 
14 -'-3.36 


17. 


1± 1.5 



Table 2: Best fit values of n and Q (quadrupolc) 
from this and previous works; the best-fit values 
of Q for n = 1 are also given. In the first column 
we give the initial of the first author (G for G96, H 
for Hinshaw et al. 1996, B for BW and T for the 
present work). The number after G indicates the 
width of the cut around the galactic equator and, 
in general, notation is the same as in G96. Re- 
sults by other authors are obtained by operating 
the so-called customary cut suggested by Bennct 
ct al. (1996). The letter w states that results were 
obtained through a weighted average of the 6 maps 
(3 frequencies, 2 channels); w2 indicates that only 
4 maps were used, neglecting 31.5 GHz outputs; 
I indicates that a linear combination of results at 
various frequencies was done, taking coefficients 
able to cancel the contribute from free-free galac- 
tic emission; — g indicates that a correction for 
foregrounds was done; results obtained consider- 
ing separately results at 3 frequencies are marked 
with * (961+961+961 component vectors). 





a 


b 


c 


S 


^0 


8.34 





-8.4 


1.89 


Pi 


18.35 


4.1 


-19.5 


1.14 


P2 


-1.574 


4.58 





2.30 


P3 


0.74 


7.4 


-8.4 


2.54 


Pi 


4.66 


-22.3 


20 


2.96 


p 5 


2.383 


-1.19 





2.64 


n° 


1.06 


-0.09 


0.38 


1.9 




6.312 


-0.1 


-0.2 


2 



Table 3: Best-fitting values of parameters in eq. 
(23) 



"custom-cut", described by Bennet et al. (1996), 
are removed. 

The level of consistency between this work and 
BW can be appreciated in Fig. 5, where the likeli- 
hood distributions on the spectral index n, found 
by us and BW, are compared. As far as the distri- 
bution on n is concerned, we reproduce even minor 
features visible in previous outputs. On the con- 
trary, our aio values tend to be slightly smaller and 
the discrepancy from BW amounts to ~ 0.64^if . 

It may be also interesting to compare our 1— 
2-3 a likelihood contours with likelihood contours 
obtainable from BW. It ought to be noticed that 
working out such contours from BW is far from 
trivial, as one has to translate their D' and D" 
parameters into n and f^A values and this requires 
a significant numerical effort. Furthermore, the 
peak likelihood curve, at various n values, ob- 
tained by BW, shows features indicating a non- 
Gaussian behaviour in respect to n. In the aio 
direction, instead, we have just the information 
that 1-cr errors (casual) are ~ 7% and only can 
we assume a Gaussian behaviour. 

However, using such information, it is possi- 
ble to work out l-2-3cr likelihood contours. In 
Fig. 6 we compare our and BW curves (thick and 
dashed lines, respectively), after displacing BW 
along the aio axis, so to have coincident top like- 
lihood points. Those interested in likelihoods ex- 
trapolated from BW results, for the case = 0, 
h = 0.5, can work them out from Fig. 6, shifting 
the aio axis by 0.64 /j,K. 

The main interest of this figure, however, con- 
cerns the actual distribution of the likelihood, once 
the non-Gaussian behaviour is fully taken into ac- 
count. Fig. 6 confirms that the distributions along 
n, in this and BW analyses, are pretty similar. 
The distributions along aio, instead, are different, 
as is expected. Substancial discrepancies already 
exist at 1-cr level and they are asymmetrical in the 
two directions of the aio axis. 

4. PS predictions and a numerical algo- 
rithm 

In the preparation of this work, as previously 
outlined, we have extensively used the public pro- 
gram CMBFAST. Besides of the transfer function 
T(k), this program provides the angular spectra 
Ci (for temperature fluctuations and polarization) , 
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normalizing all results and providing the model 
likelihood on the basis of BW relations. As previ- 
ously outlined, this implies that, for each model, 
the normalization corresponding to the best-fit to 
COBE-DMR 4-year data is selected; the model 
likelihood provided by CMBFAST is the one which 
corresponds to such best fit. 

Making use of the transfer function, a fair deal 
of large scale observables can be predicted. For 
instance, using the Press & Schechter (PS) ap- 
proach, one can evaluate the expected number 
density of galaxy clusters. As previously outlined, 
this is a fair example of why one needs to go be- 
yond the best-fit normalization. 

Using the fitting formulae worked out in the 
previous section, we have therefore built a rou- 
tine, which may be appended to CMBFAST, to 
replace to BW single normalization with possi- 
ble normalizations for top likelihood and at the 
1-2-3 a limits (this routine is available on SPOrt 
web-page, http://sport.tesre.bo.cnr.it/). Likeli- 
hood estimates are also suitably provided to com- 
pare different models, even with different normal- 
izations. 

5. Conclusions 

The main results of this work are: (i) The ex- 
pression (22) of the 2-dimensional interpolating 
"curve" yielding the likelihood distribution, for 
models with given A, when varying a\o and n. It 
is important to stress that, although the likelihood 
distribution is clearly non-Gaussian, all deviations 
from a Gaussian behaviour are fully under control 
when cubic terms are added, (ii) The expression 
(22) is then implemented by the interpolating ex- 
pressions of the Pi coefficients, given by cq. (23); 
here again, besides of outlining its practical use 
in association with Table 3 parameter values, the 
reader's attention is to be attracted to its simplic- 
ity. Notice that, given the values of the 4 param- 
eters a, b, c and S, our expressions provide the 
likelihood distribution for any critical A model in 
the a\o,n plane. That this likelihood is readily 
obtainable, for such a wide range of cases, just as- 
signing a few numerical values, is one of the results 
of this work. 

These results can be completed by an expres- 
sion allowing to compare different A models, in 
the point of the aio, n plane where they reach top 



likelihood. In Fig. 7 the top likelihood of models 
is plotted against Ha- The best-fit curve shown in 
Fig. 7 has equation: 

C = a + bni + cQf (24) 

with a = 0.98, b = 0.34, c = -1.24, S = 1.42. This 
expression is quite similar to eq. (23), yielding the 
Q,A dependence of the coefficients Pj. 

As above outlined, the class of models discussed 
here was previously treated also by BW, whose re- 
sults also concerned models with non-critical den- 
sity. If we consider only critical ACDM models, 
seemingly favoured by current observations, the 
likelihood distribution in respect to the spectral 
index n, found in BW, is the same as the one found 
in this work. In respect to normalization, instead, 
we find top likelihood for values smaller by ~ 7 % 
(see fig. 6), in respect to BW. By itself this is not 
surprising, owing to various differences between 
ours and BW analyses: first of all, the different 
techniques adopted to subtract the galactic equa- 
torial band. Such percentage can be also consid- 
ered the typical discrepancy among the amplitudes 
that different authors obtain from COBE-DMR 
data analysis, and the size of the 1-a errorbar for 
the amplitude. 

Such size of errorbar had already been sug- 
gested by BW, although no justification of their 
finding is reported in their article. Besides of a 
detailed treatment of error analysis, in this work 
we provide errors up to 3-<r, in the 3-dimensional 
space spanned by amplitude, spectral index and A, 
fully accounting for their deviation from a Gaus- 
sian behaviour. 

The above uncertainty of 7% makes perhaps 
redundant the information on non-Gaussian be- 
haviour at the 1-a level, which implies corrections 
smaller than 7 %. However, if a likelihood distribu- 
tion above 1-a is needed, the results of this work 
are to be applied; furthermore, above 1-a, devi- 
ation from Gaussian behaviour approach errorbar 
size. A typical case when likelihood distributions 
above 1-a are needed is when COBE likelihood is 
to be considered together with the likelihood dis- 
tribution worked out from other experiments. In 
this case, the joint likelihood might well reach its 
top for a model discrepant more than 1-a from 
COBE maximum likelihood, but, even if not so, 
distributions above 1-a are to be considered to 
perform a complete analysis. 
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As above outlined, our results hold almost in- 
dependently of the value of the Hubble parameter 
H. Its impact on Ci up to I = 30 is known to 
be quite limited, and we tested that, in the range 
50-70 km/s/Mpc, we do not need to specify its 
value. The same can be said for the baryonic con- 
tent of the model. Of course, both H and f2& affect 
Ci at greater I values. We have not checked our 
output against a variation of the nature of dark 
matter. If a significant massive neutrino contribu- 
tion to dark matter is present (density parameter 
Oh), some changes of C\ are expected for / ~ 200. 
Such changes are however rather small (Dodelson 
ct al. 1996), unless neutrinos have a very late 
derelativization. The situation might be differ- 
ent if hot-dark-matter with non-thermal distribu- 
tion is considered (volatile models; see Pierpaoli & 
Bonometto, 1999). In any case, however, the im- 
pact on Ci of a hot (or volatile) component, down 
to I ~ 30, is expected to be even smaller than 
the one arising from shifts. Hence, our fitting 
relations can be safely used for any value of H, 
fib and for any reasonable value of Qh originating 
from massive neutrinos. 
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A. Appendix A 



The sky temperature measured in a direction n= 6, <fi is a convolution of the temperatures around (n) 
with the beam pattern. Temperatures are then averaged inside each pixel area. Let then w plx (9, </>) be the 
filter function, that models the pixel shape. Assuming it to be symmetric around the polar axis (that we 
choose in the 9 direction), w ptx depends only on the angular coordinate and can be expanded in Legendre 
polynomials; the coefficients of the expansion then read: 



pix _ ^dSsmew^{e)p l {cose) 
w i — ■ 



/* d9 sin 9w plx {9) 



(Al) 



For COBE pixelization w plx is a circular top-hat window with angular area Q p i X = 47r/6144; let then 9 be 
such that 



d9 = 0, 



then the eq (Al) becomes 



(A2) 



w 



■ptx 



2tt 

f2 tWtt; 



pcosO 



dx Pi{x) = 



2tt 

Qrtix 



cosS , i d i (x 2 - iy 

dx 



2ir 1 r c 



dx 



dx 



ji-i 



(x 2 - iy 



dx 1 



2tt 1 
o . opn 

^ L plX ^ 1 • 



2H\ dx 1 

S-\ x 2 -iy 



dx 1 



The values of wf tx are reported in the table below. 



(A3) 



B. Appendix B 

By integrating the likelihood C over the first 4 components of the vector c, we can remove the contribu- 
tion of the monopole and dipole terms. Such integration can be easily performed in the case of Gaussian 
distributions. 

In fact, let us rewrite the expression of 

£<c »^7(lOT oxph(1/2)cTM " c1, (B1) 

using a greek index for the components from 1 to 4, and x or y for the components from 5 to 961. The 
integral of eq. (Bl) reads then 

J dc^ £(c„) = y=^== J d Cfl exp (-^c^M^ 1 ^ - ^c x M^c y - CfJt M-£c x ). (B2) 



I 


ptx 


I 




I 


pix 

w p 


I 


pix 

wf 


I 


ptx 


I 


ptx 


1 


0.9998 


6 


0.9966 


11 


0.9893 


16 


0.9780 


21 


0.9629 


26 


0.9439 


2 


0.9995 


7 


0.9954 


12 


0.9873 


17 


0.9753 


22 


0.9594 


27 


0.9397 


3 


0.9990 


8 


0.9941 


13 


0.9853 


18 


0.9724 


23 


0.9557 


28 


0.9353 


4 


0.9984 


9 


0.9927 


14 


0.9830 


19 


0.9694 


24 


0.9519 


29 


0.9308 


5 


0.9976 


10 


0.9911 


15 


0.9806 


20 


0.9662 


25 


0.9480 


30 


0.9262 



Table 4: COBE pixalization function values; see text. 
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Let be J M = M^c x , then 



/ 



x J dCf, cxp {~c tt M tl ^c-2c tt J tt ). (B3) 

Let then be = c M + 6 M(T , where 6^ is a generic matrix (4 x 4) independent from c M ; using as integration 
variables, the exponent in the integration becomes 

c^M^c + 2c M J^ = (c M - b^JvM-^Cv - b VT J T )+ 

+2(c M - J M = c^M^Cv - b^J a M~ I }c v + 
-c^M^b VT J T + b^ a J a M-^b UT J T + 2c^Jf, - 2J il b il(J J a . (B4) 
If b is the inverse of the (4 x 4) matrix (M -1 )^, the eq. (B) simplifies into 

CfiM^jj Cjy -\- 2c^J^ c^M ^ v c v J uC v J^Cf! J ^jb J 

+2c^ - 2J jl b jla J a = CpM'^Cv - J ll b iia J a . (B5) 
Henceforth, performing the integration, we have 

/ C = ^(2Tr) N deM eX ^~^ CxM * yl ° y + J ^ ffJ< ^ j dc ^ ex ^\ c ^ M ^ c ^ = 

:(2?r)2 firi^ c ^\ {CxM -y Cy + c * M A° M °v c yy = 



y/(2ir) N detM ^/(detb- 1 ) 2 



1 , 1 



(2n) N - 4 detM 



exp[--(c x M- y 1 c y )] ; (B6) 



the new covariance matrix M has ((N — 4) x (N — 4)) dimension. 

Hence, the integration of the first four components of the vector c is equivalent to replace the matrix M 
with M xy = M-y 1 + M-Jb^M-y 1 (and N with N - 4) in eq (Bl). 
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1\ 



Fig. 1. — Likelihood function for fl\ = and h = 0.5; a w is expressed in jiK. 
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Fig. 2. — Examples of comparison between the 1-2-3 a confidence levels and likelihood peaks directly 
obtained from the likelihood function (continuous lines, filled box) and the contours obtained using the 
fitting formula, eq. (22) (dashed lines, empty triangles). 




Fig. 3. — Examples of the fits between n°, a° and the parameters Pi and the analytic expression 23. 
Different symbols are related to the values of h. Discrepancies between parameters and fitting expression 
are < 1%. 
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Fig. 4. — The 1-2-3 a confidence levels on Q-n plane from 4-years COBE data, using a Sachs & Wolfe 
spectrum. The empty triangle is the top likelihood point we obtain. For the sake of comparison, the filled 
box is the top likelihood point obtained in G20+* (see text). 
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Fig. 5. Likelihood distribution against n from BW (dashed line) and our work (continuous line). 
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Fig. 6. — Confidence contours are shown for a model with f^A = 0. and h — 0.5 (thick lines) and compared 
with confidence contours obtainable from BW (dashed lines), displaced to overlap top likelihood points. This 
figure exhibits the effect of taking suitably into account the deviations from a Gaussian behaviour. 
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Fig. 7. — Top values of the likelihood for models with different The values of aio and n giving such 
likelihood are obtained from expression 24. 
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